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PROCEDE POUR L 1 ^TlKSTTON 
DES CARACTERISTIQUES D'UNE PRECIPITATION 

La presente invention concerne un procede pour 
1' estimation des caracteristiques d ! une precipitation, et en 
particulier du taux precipitant pour une precipitation solide. 

On connait dans l'Stat de la technique 1 1 utilisation des 
caracteristiques d 1 une image radar pour estiraer les 
caracteristiques pluviometriques . En particulier, le brevet 
europeen EP1049944 d§crit une technique pour 1 1 estimation de 
pluie avec un radar. Selon ce brevet, on realise les etapes 
suivantes : 

-on mesure, au moyen dudit radar bipolaire, sur un 
intervalle [ro, ri] donne de rayon de trajet par 
rapport audit radar, la phase dif f erentielle (Odp et 
la reflectivity apparente Z. selon au moins I 1 une des 
polarisations H ou V ; 

-on determine une estimation de la valeur No* 
representative de la distribution dimensionnelle des 
gouttes de pluie, a partir de la difference de la 
phase dif f erentielle entre ro et ri et a partir d'une 
integrale d'une fonction de la reflect ivite apparente 
Z, le long de 1 1 intervalle [ro, ri] ; 

-on deduit la valeur du taux de precipitation en un 
point a partir de No* et de la r6flectivite apparente 
en ce point . 

On connait egalement un proced6 pour 1' estimation de 
pluie decrit dans le brevet PCT WO03007016 dScrivant un 
procede pour 1' estimation d'un taux precipitant au moyen d'un 
radar bipolaire, caracterisS par les diff£rentes Stapes 
suivantes : 
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- on mesure au moyen dudit radar bipolaire, sur un 
intervalle [rl, rO] donne de rayon r de trajet par 
rapport audit radar, la phase dif f erentielle Odp et 
la reflectivite attenuee Z selon au moins l'une des 

5 polarisations H ou V ; 

- on determine une estimation de la valeur K(rO) de 
1 1 attenuation specif ique en rO a partir du profil de 
reflectivite attenuee ainsi mesure, ainsi qu'a partir 
de la difference de la phase dif f erentielle entre rO 

10 et rl; on determine une estimation K(r) de 

1 1 attenuation specif ique en r en fonction de 
1 1 attenuation K(r0) ainsi determin£e et du profil de 
reflectivite attenuee Z(r); on determine le taux de 
precipitant R(r) connaissant K(r) . 
15 Ges differentes solutions permettent de caracteriser des 

precipitations liquides et d'estimer de fagon precise le taux 
de pluie (en mm/h) , mais pas d'estimer les caracteristiques 
des precipitations solides telles que la neige. 

L ! invention concerne un nouveau procede visant a remedier 
20 a cet inconvenient. 

A cet effet, l 1 invention concerne selon son acception la 
plus genSrale un procede pour 1' estimation des 
caracteristiques d'une precipitation, et en particulier du 
taux precipitant pour une precipitation solide comportant une 
25 etape d 1 acquisition d'une image radar comportant au moins un 
plan vertical d f une zone de precipitation et d'un traitement 
d'un profil vertical pour fournir des signaux numeriques 
representatif s de la reflectivite selon la direction verticale 
h caracterise en ce que l'on realise une etape d' integration 
30 desdits signaux representatif s de la reflectivite pour 
d£livrer un signal representatif du profil dans le plan 
vertical du diametre moyen des particules pondere par la masse 
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de chaque particule, et une etape de determination de la 
concentration des particules solides a partir des signaux 
calcules dans les etapes precedentes. 

De preference, l'etape d' integration consiste a 
determiner la variable Z(h) 1' observable radar en mm 6 /m 3 en 
fonction de 1' altitude h k partir de ladite image radar, et a 
determiner ledit diametre moyen des particules Dm(h) par 
resolution de 1' equation : 

i^-O25k eff aDtnQr^^0 m ( 2 ) 

ou : 

■ Z est 1' observable radar a inverser en mm 6 m" 3 ; 

■ D m est en m ; 

■ a et b sont des coefficients specif iques des 
particules de type « agregats ». Le coefficient a est par 
exemple egal a 35184 et le coefficient b est egal a 3.16. 

■ k eff est le coefficient d'efficacite du processus 
d'agregation a ajuster, ce coefficient k ef£ etant par 
exemple egal ^ 0.3 

L 1 integration de (2) necessite une condition a la limite 
d 1 integration. Avantageusement , ladite condition a la limite 
d' integration est determinee afin que la valeur D ni (h) au sommet 
du nuage corresponde a la valeur predetermine du nombre total 
de particules au sommet du nuage. 

Selon un mode de realisation avantageux, on determine le 
profil du nombre total de particules n t (h) par 1' equation 
suivante : 

n T (h)=x.Z(h)/D m (h) 6 

ou x est egal a 25,4 10" 18 

Selon un autre mode de realisation, on determine le 
param£tre meteorologique N 0 (h) par 1' equation suivante : 
N 0 (h)=y.Z(h)/D ra (h) 7 
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Ou y est egal £ 102 10" 18 

Selon un troisieme mode de mise en oeuvre, on determine 
le parametre meteorologique correspondant au profil du contenu 
en glace IWC(h) [en g/m 3 ] par 1' equation suivante : 
5 IWC(h)=wZ(h)/D m (h) 3 

ou w est egal a 1,25 10~ 12 

Selon un quatrieme mode de mise en oeuvre, on determine 
le parametre meteorologique correspondant au profil du taux 
precipitant R{h) solide (mm/h equivalent fondu) par 1' equation 
10 suivante : 

R(h)=r.Z(h)/D m (h) 2 ' 35 

Ou r est egal a 4,698 10" 10 

L ! invention sera mieux comprise a la lecture de la 
description qui suit, se referant a un exemple non limitatif 
15 de realisation. 

- La figure 1 represente un exemple de profil 
vertical de Z a inverser (dans cet exemple, l'isotherme 0°C 
est au niveau du sol) . 

- La figure 2 represente la comparaison du profil 
20 de Dm resultant de 1' inversion de Z par le modele d'agregation 

avec l'estimateur classique. 

- La figure 3 represente la comparaison des 
profils de NO et nT resultant de 1' inversion de Z par le 
modele d'agregation avec l'hypothese classique et des 

25 observations. 

- La figure 4 represente le profil vertical de Z a 

inverser 

- La figure 5 represente la comparaison du profil 
de Dm resultant de 1' inversion de Z par le modele d'agregation 

30 avec l'estimateur classique. 

- La figure 6 represente la comparaison des 
profils de NO et nT resultant de 1' inversion de Z par le 
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modele d' agregation avec l'hypoth£se classique et des 
observations . 

- La figure 7 represente la comparaison du profil 
de R resultant de 1' inversion de Z par le modele d' agregation 

5 avec l'estimateur classique. 

- La figure 8 represente la sensibilite de la 
restitution du taux precipitant equivalent (mm/h) a la loi de 
densite des particules p(D) °c D-y. 

- La figure 9 represente la sensibilite de la 
10 restitution de la vitesse Doppler verticale a la loi de 

densite des particules p(D) « D-y. 

Le procede s ' applique essentiellement aux precipitations 
stratiformes . II considere que les noyaux gla<?ogenes ne sont 
actives qu'i. temperature fortement negative, c'est-a-dire au 
15 sommet du nuage. Les cristaux de glace formes 3. haute altitude 
sedimentent et grossissent au cours de leur chute soit par 
sublimation de la vapeur d'eau saturante ambiante, soit par 
collection et givrage de gouttelettes d'eau nuageuse 
surfondues, soit par agregation au hasard de leurs collisions 
20 avec d'autres cristaux de glace. Des trois processus de 
croissance, seule 1' agregation change la concentration en 
particules de glace. Le procede de caracterisation repose 
essentiellement sur une description simplifiee du mecanisme 
d'agregation. L" invention est basee sur une classe de procedes 
25 « profileurs », en ce sens qu'elle inverse le profil vertical 
de reflectivity mesure dans la glace, pour en tirer le profil 
vertical du taux precipitant solide. 

Les etapes de la methode d' inversion sont les suivantes : 
1- La distribution dimensionnelle des particules 
30 exprim^e en « diametre equivalent fondu », est supposee 

exponentielle i.e. : 

N(D)=N 0 exp(-4D/Dm) (1) 
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ou N(D) est la concentration en particules par m 3 et 
par intervalle de diamStre, et N 0 et D m sont les deux 
parametres qui caracterisent la distribution. 

2- On determine le sommet et la base h min de la 
couche de precipitation solide. 

a. h max est 1' altitude maximum du profil de 
reflectivite mesurg Z{h) . 

b . h mia est soit 1' altitude de l'isotherme 0°C si 
la temperature au sol est positive, soit le niveau du 
sol si la temperature au sol est negative. 

3- On determine alors le profil du paramdtre D m (h) 
entre h max et h min en resolvant 1' equation 
dif f erentielle : 

i^^025k eff aDtn0^Z+^^)) m ( 2 ) 

ou : 

Z est I 7 observable radar a inverser en 
mm 6 m~ 3 ; 

D m est en m ; 

a et Jb sont des coefficients specifiques 
des particules de type « agrggats », egaux 
respectivement a 35184 et 3.16 & partir des 
observations de Locatelli et Hobbs (1974) ; 

k eff est le coefficient d'efficacite du 
processus d'agr£gation a ajuster (la valeur k eff 
= 0.3 semble correcte) . 

4- L' integration de (2) s'effectue £ partir du 
haut, ou la condition & la limite s'exprime en fixant le 
nombre total de particules n T (ou nombre de noyaux 
glagogenes actives au sommet de nuage) . On peut prendre 
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n r (h max )= 10 6 trf 3 , ce qui permet d'exprimer la condition a 
la limite D^h^) conime : 

Dm(H ma *)=25,4 Iff 18 (Z(h max /m{h m(a )) 1/6 (3 ) 

5- Le profil D m (h) entre h max et h min une fois 
5 determine, on calcule les profils des autres parametres 

d'interet par les expressions suivantes : 

a. Profil de N 0 : N 0 (h) =102 . 10" 12 Z (h) /D B (h) 7 

b. Profil du nombre total de particules n T (h) 
[en nf 3 ] : N T (h) =102 . 10" 12 Z (h) /D m (h) 6 

10 c. Profil du contenu en glace IWC(h) [en g/m 3 ] 

: IWC(h)=l / 25.10" 12 Z(h)/D m (h) 3 

d. Profil du taux precipitant R(h) solide 
(mm/h equivalent fondu) 

En utilisant la loi de vitesse terminale de 
15 chute determin£e par Locatelli and Hobbs pour des 

agregats: [v r =107.6 D 0 * 65 (D in m) ] # R(h) s'exprime 
par : 

R (h) =4 , 698 . 10" 10 Z (h) /D ra (h) 2 ' 35 
La description qui suit concerne une autre version du 
20 procede de traitement pour la determination du taux 
precipitant . 

Pour la pluie, 1 ' algori thme utilise est 
l'algorithme ZPHI objet des brevets FR9800714 et FR0109206. 

Pour les precipitations solides, le nouvel 
25 algori thme qui vient d'etre developpe est 1' objet de la 
description qui suit. 

Cet algorithme d' estimation du taux precipitant 
pour une precipitation solide fait partie de la classe des 
algorithmes « profileurs » t en ce sens qu'il inverse le profil 
30 vertical de reflect ivite mesure dans la glace, pour en tirer 
le profil vertical du taux precipitant solide. 
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L ' algorithme s' applique essentiellement aux 
precipitations stratif ormes . II considdre que les noyaux 
glagogenes ne sont actives qu'a temperature fortement 
negative, c'est-a-dire au sommet du nuage. Les cristaux de 
glace formes a haute altitude sedimentent et grossissent au 
cours de leur chute soit par sublimation de la vapeur d'eau 
saturante ambiante, soit par collection et givrage de 
gouttelettes d'eau nuageuse surfondues, soit par agregation au 
hasard de leurs collisions avec d'autres cristaux de glace. 
Des trois processus de croissance, seule 1' agregation change 
la concentration en particules de glace. L' algorithme repose 
essentiellement sur une description simplifiee du mecanisme 
d'agregation. Les etapes de la m<§thode d' inversion sont les 
suivantes : 

La distribution dimensionnelle des particules 
exprimee en « diametre equivalent fondu », est supposee 
exponentielle i.e. : 

N(D)=Noexp(-4D/n») , (1) 
ou: 

• D est le diametre equivalent fondu de la particule 
de glace, 

• N(D) est la concentration en particules par m 3 et 
par intervalle de diametre, 

• et N 0 et D m sont les deux parametres qui 
caracterisent la distribution. 

On determine le sommet h^ et la base ^ de la 
couche de precipitation solide. 

• h max est l 1 altitude maximale du profil de 
ref lectivite mesure Z (h) . 
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• est soit 1' altitude de 1 ' isotherme 0°C si 

la temperature au sol est positive, soit le 
niveau du sol si la temperature au sol est 
negative . 

On determine alors le profil du parametre D m (h) 
entre et en resolvant 1' Equation dif f erentielle , dont 

la solution peut etre determinee analytiquement : 



d Dm _ A d KffW)" t-ia-2) 1 az 

"aT" r(4W) d(a-l) LDm {a-\)Zdh Um 



(2) 



10 

ou: 



J(b 9 d)= rexp(-4x)x b \*\c d -x' d ]xp(-4x i )dx'dx decrit la 
Jo Jo 

frequence de collision 
15 • Z est 1' observable radar a inverser en mm 6 m" 3 ; 

• D m est en m ; 

• a, b, c et d sont des coefficients dependant de la loi 
de densite F (D) des particules de glace, supposee 

20 varier en D~ r (ou 7 peut prendre des valeurs entre 0.25 

et 1.1 suivant le type des particules). Ces 
coefficients sont tires de la theorie de Mitchell 
(Journal of Atmospheric Sciences, 53, 12, 1996) pour 
representer les lois puissance suivantes : 



25 



v t (D) = cD d (ou v t est la vitesse terminale de 



chute de la particule de diametre equivalent fondu D) 

A = aD b (oil A est la section ef f icace de la 
particule de diametre equivalent fondu D) 

30 
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Avec des coefficients dependant de la frequence et 
d'une fonction de (D) , repr£sentant par une loi puissance la 
relation entre la reflectivity (parametre mesure) , la 
concentration des particules et leur diametre equivalent 
moyen : 

Z = dN 0 D m d 

k eff est le coefficient d'efficacite du processus 
d'agregation a ajuster (la valeur k e£f = 0.3 semble correcte 
juste au-dessus de l'isotherme 0°C) . 

L' integration de (2) s'effectue a partir du haut, 
ou la condition a la liraite s'exprime en fixant le nombre 
total de particules n T (ou nombre de noyaux glagogenes actives 
au sommet de nuage) . On peut prendre n T (h^)^ 10 6 m~ 3 , ce qui 
permet d'exprimer la condition a la limite D m (h max ) comme : 



1 Z(*max) 



l 



(3) 



Le profil D m (h) entre h max et h min une fois 
determine, on calcule les profils des autres parametres 
d'int§ret par les expressions suivantes : 

Profil de N 0 : 



N 0 {hy 



dD m (h) d 



Profil du nombre total de particules n T (h) [en m" 3 ] 
n T (h) = 025N 0 (h)D m (.h) 

Profil du contenu en glace IWC{h) [en g/m 3 ] 
IWC(K) = (1.2272)10 4 N 0 (h)D m (h) 4 
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Profil du taux precipitant R(h) solide (mm/h 
equivalent fondu) : 

R{h) = (1.885)10 6 C ^f N 0 (h)D m (h)^ 

5 Profil vertical de la vitesse Doppler V D (h) (en 

m/s) , donnee par une loi puissance dependant de la densite 
_(D) : 

V D (h) = pD m {h)* 

10 

Du fait de la dependance des coefficients vis-a- 
vis de la loi de densite F (D) , la restitution du taux 
precipitant depend de fagon critique du parametre y, comme en 
temoigne la figure 7. Cette loi de densite constitue done un 

15 paramdtre clef. 

En operant en tir vertical, on peut mesurer la 
vitesse Doppler V D (h), et la comparer au V D (h) calcule par 
l'algorithme. Le profil de vitesse Doppler est discriminant, 
comme en temoigne la figure 8. Par ajustements successifs, on 

20 peut ainsi determiner le parametre y de la loi F(D) pour 
lequel la coherence entre V D (h) calculee et V D (h) mesuree est 
realisee. La figure 4 represente un exemple de profil inverse 
par l'algorithme (dans cet exemple, I'isotherme 0°C est au 
niveau du sol) . La figure 5 represente les courbes resultant 

25 de 1' application du moddle d'agregation et du modele 
classique, permettant la comparaison du profil de D m resultant 
de 1' inversion de Z par le modele d'agregation avec 
l'estimateur classique. 

La figure 6 represente la courbe des profils en 

30 fonction de 1' altitude (m) au-dessus de I'isotherme 0°C. Elle 
conduit a la comparaison des profils de N 0 et n T resultant de 
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1' inversion de Z par le modele d'agregation avec l'hypothese 
classique et des observations. La figure 7 represente la 
comparaison du profil de R resultant de 1' inversion de Z par 
le modele d'agregation avec l'estimateur classique. La figure 
5 8 represente la sensibilite de la restitution du taux 
precipitant equivalent (mm/h) a la loi de densite des 
particules p(D) oc D' y . La figure 9 represente la sensibilite de 
la restitution de la vitesse Doppler verticale a la loi de 
densite des particules p (D) oc D" Y . 

10 



WO 2005/040853 



13 

RE VEND I CAT I ONS 
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1 - Procede pour 1 1 estimation des caracteristiques d'une 
precipitation, et en particulier du taux precipitant pour une 
precipitation solide comportant une etape d' acquisition d'une 
image radar comportant au moins un plan vertical d'une zone de 
precipitation et d'un traitement d'un profil vertical pour 
fournir des signaux numeriques representatif s de la 
reflectivite selon la direction verticale z caract£rise en ce 
que l'on realise une etape d' integration desdits signaux 
representatif s de la reflectivite pour delivrer un signal 
representatif du profil dans le plan vertical du diametre 
moyen des particules ponder^ par la masse de chaque particule, 
et une etape de determination de la concentration des 
particules solides a partir des signaux calcules dans les 
etapes precedentes . 

2 - Procede selon la revendication 1, caracterise en ce 
que 1' etape d' integration consiste a determiner la variable 
Z(h) 1' observable radar en mm 6 /m 3 en fonction de 1' altitude h a 
partir de ladite image radar, et a determiner ledit diametre 
moyen des particules Dm(h) par resolution de 1' equation : 

^-025A: c# aD^0-^Z+|^||)) ;w (2 ) 

ou : 

■ Z est 1' observable radar a inverser en 
mm 6 nf 3 ; 

■ D ra est en m ; 

■ a et b sont des coefficients 
specifiques des particules de type 
« agregats » . 
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■ k eff est le coefficient d'efficacite 
du processus d'agregation a ajuster. 

3 - Procede selon la revendication 2, caracterise en ce 
5 que le coefficient k eff est £gal a 0.3. 

4 - Procede selon la revendication 2 ou 3, caracterise en 
ce que le coefficient a est egal a 35184. 

10 5 - Procede selon la revendication 2, 3 ou 4, caracterise 

en ce que le coefficient b est egal & 3.16. 

6 - Procede selon l'une quelconque des revendi cat ions 2 a 
4, caracterise en ce que ladite constante d' integration est 

15 determinee afin que la valeur D m (h) au sommet du nuage 
corresponde a ia valeur predeterminee du nombre total de 
particules au sommet du nuage. 

7 - Procede selon l'une quelconque des revendications 
20 precedentes, caracteris6 en ce que l'on determine le profil du 

nombre total de particules n t (h) par 1' equation suivante : 
n T (h)=x.Z(h)/D m (h) 6 . 

8 - Proced£ selon la revendication precedente, 
25 caracterise en ce que x est egal a 25,4 10~ 1{ \ 

9 „ Proced§ selon l'une quelconque des revendications 
precedentes, caracterise en ce que 1'on determine le parametre 
meteorologique N 0 (h) par liquation suivante : 

30 N 0 (h)=y.Z(h)/D m (h) 7 . 
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10 - Procede selon la revendication prScedente, 
caracterise en ce que y est egal a 102 10' 18 . 

11 - Procede selon l'une quelconque des revendications 
precedentes, caracterise en ce que l'on determine le parametre 
meteorologique correspondant au profil du contenu en glace 
IWC(h) [en g/m 3 ] par 1' equation suivante : 

IWC(h)=wZ(h)/D n (h) 3 . 

12 - Procede selon la revendication precedente, 
caracterise en ce que w est egal a 1,25 10" 12 . 

13 - Procede selon l'une quelconque des revendications 
precedentes, caracterise en ce que l'on determine le parametre 
meteorologique correspondant au profil du taux precipitant 
R(h) solide (mm/h equivalent fondu) par 1' equation suivante : 

R(h)=r.Z(h) /D m (h) 2 ' 35 . 

14 - Procede selon la revendication precedente, 
caracterise en ce que r est egal k 4,698 10" 10 . 
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Fig 4 
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Fig 6 
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Fig. 8 
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